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I ■ ABSTRACT 

O ; 

[ Aims. We combine measurements of weak gravitational lensing from the CFHTLS-Wide survey, supernovae la from CFHT SNLS and 

CMB anisotropies from WMAP5 to obtain joint constraints on cosmological parameters, in particular, the dark-energy equation-of- 
state parameter w. We assess the influence of systematics in the data on the results and look for possible correlations with cosmological 
parameters. 

, Methods. We implemented an MCMC algorithm to sample the parameter space of a flat CDM model with a dark-energy component 

of constant w. Systematics in the data are parametrised and included in the analysis. We determine the influence of photometric 

1 calibration of SNIa data on cosmological results by calculating the response of the distance modulus to photometric zero-point 

^vq ' variations. The weak lensing data set is tested for anomalous field-to-field variations and a systematic shape measurement bias for 

high-redshift galaxies. 

. Results. Ignoring photometric uncertainties for SNLS biases cosmological parameters by at most 20% of the statistical errors, using 

supernovae alone; the parameter uncertainties are underestimated by 10%. The weak-lensing field-to-field variance between 1 deg 2 - 
MegaCam pointings is 5%-15% higher than predicted from TV-body simulations. We find no bias in the lensing signal at high redshift, 
within the framework of a simple model, and marginalising over cosmological parameters. Assuming a systematic underestimation 
of the lensing signal, the normalisation cr 8 increases by up to 8%. Combining all three probes we obtain -0.10 < 1 + w < 0.06 at 68% 
confidence (-0.18 < 1 + w < 0. 12 at 95%), including systematic errors. Our results are therefore consistent with the cosmological 
!L* constant A. Systematics in the data increase the error bars by up to 35%; the best-fit values change by less than 0.15cr. 



Key words. Cosmology: observations - Cosmological parameters - Methods: statistical 

1 . Introduction cosmic shear data dFu et al.l f2008. hereafter F08) and supernova 

(SN) la data fro m the Supernova Legacy Survey (SNLS, 

The Canada-France-Hawaii-Telescope Legacy Survey |Astier et all 12001 hereafter A06). We complement these with 

(CFHTLS3) has yielded interesting constraints on cosmo- CMB measurem ents from the Wilkins on Microwave Anisotropy 

logical parameters using different probes, for example weak p ro be (WMAP5,|HinshawLiLaDllQQ9l hereafter H09). 
gravitational lensing and supernovae of Type la. In this paper 

we combine two recent measurements from the CFHTLS: the ™ i_ j j r ^ T c i_- i_ i_- i_ 

There are hundreds of type-la supernovae for which high- 

quality observations are available. Many surveys spanning wide 

Send offprint requests to: Martin Kilbinger, e-mail: redshift ranges and using different telescopes and search strate- 

kilbinger@iap.fr gj es can ^ e accessec [ a recent compilation of various sur- 

* Based on observations obtained with MegaPrme/MegaCam, a joint ve ys and the derived cosmological constraints can be found in 

project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii | KoW alski et all (120081) . Despite this impressive number of avail- 

Telescope (CFHT) which is operated by the National Research Council ui u- * u u . i ..u cmt c 

„„„ F r A , , „ . , ,, TT . r. , able objects we have chosen to restrict ourselves to the SNLS 

(NRC) ol Canada, the Institut National des Sciences de 1 Univers ol the , ,. . .„ „ r 

Centre National de la Recherche Scientifique (CNRS) of France, and sam P le - The resulting constraints will suffer from a greater sta- 

the University of Hawaii. This work is based in part on data products tlstlcal uncertainty; however, the systematic errors are better un- 

produced at Terapix and the Canadian Astronomy Data Centre as part derstood by using a single, homogeneous survey. By providing 

of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative imaging with a single telescope and camera for example, SNLS 

project of NRC and CNRS. allows a common photometric calibration strategy for the whole 

1 http://www.cfht.hawaii.edu/Science/CFHTLS survey. 
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The lensing data presented in F08 might suffer from sys- 
tematic errors in the shear calibration and shape measurement. 
In this work, we performed two tests for the presence of sys- 
tematics. First, we compared the variance between individual 
MegaCam pointings to simulations. Second, we quantified a po- 
tential systematic underestimation of the lensing signal at high 
redshifts. The influence on cosmological parameters was esti- 
mated. Due to the unknown origin of those systematics, the anal- 
ysis must remain less rigorous than for the SNIa case. It is by no 
means comprehensive since we only quantify the effect of those 
systematics and do not investigate their origin, which is beyond 
the scope of this paper. 

The analysis of the five-year WMAP data has yielded im- 
pressive results on many cosmological parameters from a great 
number of models dPunklev et all 120091: iHinshaw et ail 120091: 
iKomatsu et al.|[20Q9h . Thanks to the rich features in the angular 
power spectrum, many cosmological parameters can be deter- 
mined with high precision. However, degeneracies between pa- 
rameters remain, in particular, for dark-energy models. In order 
to lift these degeneracies, measurements of CMB anisotropics 
have to be complemented with other probes. 

The WMAP5 team has combined their data with other 
probes which are sensitive only to the geometry of the Universe, 
i.e. SNIa and baryonic acoustic oscillations (BAO). Other teams 
have incl uded probes of stru cture growth, like counts of X-ray 
cluste rs dMantz et all 120081) . SDSS and Lyor-forest dXia et alj 
2008). In this work, we combine WMAP5 with weak gravita- 
tional lensing which is sensitive to both the geometry and the 
growth of structure. Probing both domains will allow future sur- 
veys to distinguish between dark energy and modified gravity as 
a possible cause for the present acceleration of the Universe (e.g. 
iJain & Zhangll2 008). Such a test will be feasible with upcoming 
and proposed surveys such as KIDS0, DESQ LSSlfl JDErvfl and 
EuclicQ 

The first-year WMAP data has alr eady been combined w ith 
weak lensing, using the RCS sur vey (|Hoekstra et alj|2002l) . to 
improve constraints on Q, m and <x 8 dContaldi et alj2003 ). In F08, 
cosmic shear has been supplemented with the third-year WMAP 
data. In this work we extend the latter, simple analysis by drop- 
ping the (not well motivated) priors on some parameters for lens- 
ing (£l b , n s , h). 

The 'concordance' flat ACDM model of cosmology provides 
an excellent fit to WMAP5 and most other probes of the geome- 
try and large-scale structure of the Universe. It only contains six 
free parameters, Q.^, Q m , t, n s , h and cr% (or functions thereof). 
This model assumes a cosmological constant A as the cause 
of the observed accelerated expansion of the Universe today. 
The cosmological constant has yet evaded all plausible physi- 
cal explanations of its nature and origin and, further, brings in 
problems of fine-tuning and coincidence. A signature of some 
'dark energy' beyond a simple cosmological constant might be 
an equation-of-state (eos) parameter w = pipe 2 which deviates 
from the vacuum energy value of -1. We therefore extend the 
concordance six-parameter model by including the dark-energy 
eos parameter w. Although this vvCDM model with a constant 
w is not better motivated physically than the cosmological con- 
stant, a significant observed deviation from w = — 1 will defi- 
nitely be an indication for new physics. Moreover, the data at 



2 http://www.astro-wise.org/projects/KIDS 

3 https://www.darkenergysurvey.org 

4 www.lsst.org 

5 http://universe.nasa.gov/program/probes/jdem.html 

6 http://sci.esa.int/science-e/www/area/index. cfm?fareaid= 102 



present are not good enough to constrain more dark energy pa- 
rameters in a general way. Thi s is only feasible fo r very specific 
models, e.g. quintessence (see ISchimd et alj|2007l) or early dark 
energy (Bartelmann et al. 2006; Franci s et alJl2009h . and will be 
subject of a future work using the CFHTLS data (Tereno et al. in 
prep.). 

One goal of this paper is to focus on systematic errors, nui- 
sance parameters and their interplay with cosmological param- 
eters. For the SNIa data, apart from the usual light-curve pa- 
rameters, we take into account photometric calibration errors. A 
joint Bayesian analysis including systematic and cosmological 
parameters is done and correlations are revealed. This will be 
mandatory for future surveys with statistical errors which will 
be smaller by several orders of magnitude as compared to today. 
To decrease systematics further and further is very challenging, 
and technical limitations might set a barrier to this endeavour. 
Therefore, it is important to quantify the effect of systematics 
and nuisance factors on cosmological constraints. 

This paper is organised as follows. Section|2]contains a brief 
description of the data together with the likelihood used later in 
the analysis. It also discusses sources of systematics in the data 
and their correlations with cosmological parameters. In Sect. [3] 
we define the cosmological model tested in this work, introduce 
our implementation of the Monte Carlo Markov Chain (MCMC) 
technique and present the cosmological results of the analysis. 
We conclude with a discussion in Sect. [4] and an outlook in 
Sect.H 



2. Data and method 

2.1. CFHTLS-Wide cosmic shear 

We use the cosmic shear results from the CFHTLS-Wide 3 rd 
year data release (T0003, iFu et alj I2008I) . On 57 square de- 
grees (35 sq. deg. effective area), about 2 x 10 galaxies with 
/AB-magnitudes between 21.5 and 24.5 were imaged. The data, 
reduction analyses and shear pipeline are des cribed in detail 
in F0 8. We use the aperture-mass dispersion (Schneide r et al.l 
1998) measured between 2 and 230 arc minutes. Due to the com- 
pensated nature of its filter, this second-order measure is least 
sensitive to large-scale systmatics in the data. The source red- 
shift distribution is obtained by using the CFHTLS -Deep p(z) 
(Ilbert et al.1 12006) and by rescaling it according to the Wide i^B 
magnitue distribution and weak-lensing galaxy weights. 

To ascertain that the quality and reliability of the shear mea- 
surements are sufficient for this work, we perform further tests of 
the data beyond what has been done in F08. This is addressed in 
Sects. [2TT31 and [2~1.2| where we assess the importance of poten- 
tial systematics for the current data, and estimate their influence 
on inferred cosmological parameters. 

As in F08 the log-likelihood is In L = -^ 2 /2-ln |C|/2+const, 
where the x 2 is modelled as 

= Z (Kp (0i) L - K>(^)» i c ~% 

x (K^L-K^'^))- cd 

The predicted aperture-mass dispersion given by a model pa- 
rameter vector p is fitted to (M^ p ) b s measured at angular scales 

9j. The covariance C of (M*) b S is the one used in F08 and con- 
tains shape noise, (non-Gaussian) cosmic variance and the resid- 
ual B-mode. Those parts of the covariance which depend on the 
shear correlation (mixed and cosmic variance terms) are calcu- 
lated using a theoretical model of the large-scale structure and 
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therefore depend on cosmological parameters. We ignore this 
dependence and keep the covariance constant, corresponding to 
the fiducial cosmology as in F08. This biases the posterior confi- 
dence regions, but the effect is weak over t he region of param eter 
space permitted by CMB and lensing, see lEifler et al.l (120081) for 
a detailed discussion. Moreover, we can drop the term In \C\/2 
in the log-likelihood because its paramete r-dependence mani- 
fests i tself only for very small survey areas dKilbinger & Munshl 
120061) . 

As usual, the following relation between the aperture- 
mass dispersion and the weak lensing power spectrum holds 
dSchneider et al.ll 19981) . 
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The lensing power spectrum is a projection of the 3d matter- 
densi ty power spec trum Pg, weighted by the redshift distribution 
P (X) (lKaiseri ll992h. 
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For the non-line ar evolution of the power spectrum, the fit- 
ting formula of ISmith et all d2003l) is used. Although tested 
for ACDM models, it provides r easonable good fits to wCDM 
cosmologies as well dMal 12007b - The accuracy of any non- 
linear fitting function i s limited; when using the a nsatz of 
iPeacock & Doddsl (1 19961) instead of ISmith et all d2003l) . the re- 
sulting best-fit cr 8 differs by 2% (F08). 

We parametrise the redshift distribution using the function 



P(Z) oc 

z + c 
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p(z)dz = 1, 



(4) 



which we fit to the obtained redshift histogram. The correspond- 
ing^ 2 is 



S 



[n t - p(zd] 



(5) 



Here, n, is the normalised number of galaxies in the /* red- 
shift bin and p(zi) the fitting function, evaluated at the red- 
shift bin centre. The uncertainty cr, of «, contains Poisson noise, 
phot o-z uncertainty and c osmic variance, as described in F08 
and lBeniamin et al] d2007l) . The sum in eq. (0 extends over the 
redshift range < z < 2.5. In this range, cosmic variance is the 
dominant uncertainty. As in F08 we neglect the cross-correlation 
between different bins. 



2.1 .1 . Systematics in the lensing data 

In the following two subsections we address the question of po- 
tential residual systematics in the cosmic shear data. We estimate 
the influence of some systematics on cosmological constraints. 
The shear catalogue used here has been extensively tested in 
F08. Internal consistency checks have been performed involv- 
ing the comparison between two data reduction pipelines, resid- 
ual B-modes and cross-correlations between galaxies and stars. 
In addition, the shape measurement pipeli ne has been calibrated 
with the STEP1 and STEP2 simulations dHevmans et al.ll2006t 



Massev et al. 2007), and a bias of less than 3% has been deter- 
mined. 

Despite that, there are indications of remaining systemat- 
ics which do not manifest themselves in the merged shear cat- 
alogue. For example, there are variations of the shear signal be- 
tween individual MegaCam pointings which vanish on average. 
Moreover, with the help of photometric redshifts covering parts 
of the survey, a problem with the redshift-scaling of the shear 
correlation becomes apparent. 

2.1.2. Variance between CFHTLS pointings 

In F08 it was shown the systematics are small globally. The B- 
mode is consistent with zero on most scales (note however a 
significant detection around 100 arc minutes). The level of sys- 
tematics is low on average, e.g. the residual cross-correlation be- 
tween uncorrected star ellipticities and galaxy shapes after PSF 
correction. There might however be systematics present on in- 
dividual fields. Even if they vanish on average they are an addi- 
tional noise source and increase the measurement error bars. 

Here, we address the question whether individual pointings 
show anomalous variations which might be due to incorrect 
shape measurements or PSF correction. The field-to-field vari- 
ance of the lensing signal is compared with numerical simula- 
tions. 

We measure the shear aperture-mass dispersion on all 57 
fields, each one corresponding to a 1 deg 2 -MegaCam image, and 
calculate the field-to-field variance. W e use A^-body simulations 
dYahagi & Yoshii|[200U Yahagi 2005) through which we shoot 
light rays to obtain shear maps (Hamana & Mellier 2001) using 
a redshift distribution corresponding to the mean n(z) estimated 
for the observations. The underlying model is a flat WMAP3- 
cosmology (Q m = 0.27, Q b = 0.044, n s = 0.95, h = 0.71, cr 8 = 
0.77). 

An important issue are the correlations between CFHTLS- 
Wide pointings due to large-scale structure. Unfortunately, we 
cannot simulate the exact survey geometry of the CFHTLS-Wide 
T0003 survey which extends up to 8 degrees angular separation, 
since our simulations are only 3° x 3° in size. We therefore use 
several (independent) ray-tracing simulations to cover the ob- 
served fields; three for Wl, one for W2 and four for W3. This 
results in less-correlated fields and a smaller cosmic variance 
compared to the observations. We therefore expect the simula- 
tions to slightly underestimate the observed field-to-field vari- 
ance. The galaxies on each simulated pointing are distributed 
homogeneously, we do not take into account masking. We do 
however simulate the varying number density between fields by 
applying the observed number densities to the simulations. 

From Fig. Q]we see that the observed E-mode variance of 
(M 2 p ) is 10-20% higher than the predicted value, on scales be- 
tween 1 and 30 arc minutes. The B-mode variance is in very 
good agreement, with a tendency to be slightly lower than pre- 
dicted. From the fact that the B-mode variance is of the same 
order as the E-mode, one sees that the considered angular scales 
are shape-noise and not cosmic-variance dominated. The power- 
law shape of the variance is another indication for this. The low- 
level oscillations on scales on the order of the (M 2 p )-correlation 
length are probably due to noise. 

The observed top-hat field-to-field variance is substantially 
greater than the one from simulations, in particular at large angu- 
lar scales. This might be due to its higher sensitivity to residuals 
on large scales compared to the aperture-mass. Note however 
that because of the high correlations between different angular 
scales for top-hat, the significance of any discrepancy is hard 



4 



M. Kilbinger et al.: Dark-energy constraints from weak lensing, SNIa and CMB 



Aperture-mass dispersion 



0> o 
> c: 

132 
m <£ 
ccs 



1e-04 



1e-05 



1e-06 
0.4 

0.2 



-0.2 




Simulations E-mode 
Simulations B-mode 
Wide T0003 E-mode 
Wide T0003 B-mode 
F08 merged cat. 




6 [arcmin] 
Top-hat shear rms 



1e-04 



m 1 e-05 



1e-06 
0.6 

0.4 

0.2 






Simulations E-mode 
Simulations B-mode 
Wide T0003 E-mode 
Wide T0003 B-mode 
F08 merged cat. 




6 [arcmin] 



10 



Fig. 1. Variance of (M* p ) from the numerical simulations (solid 
lines) and CFH TLS data (dash ed). The rms used for the merged 
catalogue (from lFuet al.ll2.Q08h . which we use for the cosmolog- 
ical constraints, is the dotted line. 



to interpret. Due to the small number of simulations we cannot 
attempt an error estimate on the field-to-field variance. 

Note that the variance for the merged catalogue used in F08, 
and this work for the cosmological constraints, is roughly a fac- 
tor of four less than the one from individual fields. 

The greater observed field-to-field fluctuations could be a 
sign for uncorrected residual systematics. We have not included 
in the simulations the varying mean redshift due to the dif- 
ferent numbers of exposures for each field in the final stacks. 
The variance in limiting magnitude between fields is 0.25. With 
the empiri cal law between lim i ting m agnitude and mean red- 
shift from van Waerbek e et al. (2006) this translates into a z- 
variance of 0.03. Using the approximate relation from linear the- 
ory (Ml p ) oc z l ' 5 we find an additional expected variance of about 
4.5%. We conclude that the observed E-mode field-to-field vari- 
ance is higher than predicted by at most 5%-15%. It is difficult 
to assess the influence of this additional error on the complete 
galaxy catalogue used here. First, we do not know how scales 
larger than 30 arc minutes are affected. Second, the shear corre- 
lation used to constrain cosmology is calculated with many more 
galaxy pairs than in this field-to-field analysis, with a large num- 
ber of pairs stemming from different MegaCam pointings. 



2.1 .3. Systematics in the redshift-scaling of the shear signal 

Preliminary analysis of weak lensing tomography of the 
CFHTLS-Wide shows that beyond a redshift of about unity, the 
cosmic shear signal does not increa se as e xpected, but instead is 
systematically underestimated dFull2008l) . For this analysis we 
use the fourth CFHTLS data release (T0004) which provides 
five-band photometry on 35 square degrees. Photometric red- 
shifts for each galaxy have been obtained using the template- 
fitting code Le Phare (ICoupon et alj [2009). The quality of the 
new p hoto-z's is c ompatible to those used in F08 (T0003 data re- 
lease, Ilbert et al.l l2006). An improvment has been obtained due 
to additional and larger spectroscopic samples. We take a sub-set 
of T0004 consisting of a contiguous area of 19 square degrees. 
The measured galaxy shapes are the same as in F08. 

As an illustration, we show the top-hat shear variance mea- 
sured at 5 and 25 arc minutes for various redshift bins, see 
Table Q] The measured values, corresponding to 19 square de- 
grees, are compared to predictions using a flat ACDM model 
with Q m = 0.25 and <x 8 = 0.8. For the highest bin with redshifts 
above 0.95, the measured values are clearly inconsistent with the 
predictions. 



Table 1. Shear top-hat variance at 5' and 25' and for five redshift 
ranges. The error bars are Poisson noise. 





<lr 2 (5')|)/io 5 


<| 7 2 (25')l>/10 5 


z-range 


CFHTLS 


Prediction 


CFHTLS 


Prediction 


0.0 ... 


0.65 


-0.3 ±2.0 


0.95 


0.2 ±0.5 


0.29 


0.0 ... 


2.5 


1.9 ±0.9 


2.38 


0.7 ± 0.2 


0.73 


0.65 . . . 


0.95 


2.9 ±2.2 


4.12 


1.1 ±0.5 


1.27 


0.65 . . . 


2.5 


3.3 ± 1.4 


4.77 


0.9 ± 0.3 


1.48 


0.95 . . . 


2.5 


-1.7 ±3.7 


7.46 


0.04 ± 0.9 


2.33 



We have to assume that this anomalous redshift-scaling is 
also present in the shear catalogue used in this paper for the cos- 
mology constraints. The redshifts for this preliminary analysis 
are taken from the CFHTLS Wide and are therefore slightly less 
accurate than the T0003-Deep ones, used to infer the F08 red- 
shift distribution. However, the former do not suffer from cos- 
mic variance. The calibration with spectroscopic redshifts as- 
sures their reliability for fog < 24. The depth of the F08 shear 
galaxies is only half a magnitude higher, and those faint galaxies 
are down-weighted in the weak lensing analysis. We therefore 
assume that the problem is mainly due to the shape measurement 
and not to the photometric redshifts. One reason for this could 
be a shear calibration bias which depends on galaxy properties 
that are a function of redshift. A bias in the shape measurement 
which is a function of galaxy size, magnitude or galaxy signal- 
to-noise will affect low- and high-z galaxies differently. 

For this work we choose a simple toy model to parametrise 
this potent ial underestimation of the redshift-scaling of the shear 
signal, see Lesgourgues et al. (2007J) for a similar approach. This 
decrease of the lensing signal corresponds to an effective lens- 
ing efficiency which is lower than expected at high redshifts. To 
model this we multiply the redshift distribution n(z) with a con- 
stant co > 0, for redshifts z > Zo- We chose zo to be 1.0. en = 1 is 
the unbiased case with no degradation of the shear signal, co = 
is the pessimistic case where there is no shear signal for z > Zo- 
We allow values greater than unity corresponding to an overes- 
timation of the shear signal. Note that we use the original, unal- 
tered redshift distribution when fitting the n(z) histogram (eq.|5|. 

Using a model with fixed Qb and n s and marginalising over 
h, cr 8 and redshift parameters we measure a value of cq consistent 
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Fig. 2. Left panel: 2d marginals (68% and 95%) for £2 m and <x 8 
for the two cases of including the high-z calibration bias (solid 
lines) and ignoring it (dashed curves), respectively. Right panel: 
Id marginal likelihood for cq. The vertical lines indicate mean 
(dashed) and 68% confidence intervals (dotted). 



the PSF correction have been made in F08, and extensive studies 
going into much more detail as presented here will be published 
soon (van Waerbeke et al. in prep.). 

2.2. SNLS supernovae la 

The supernova data set corresponds to the first data release of 
the SNLS, see A06. We use 105 supernovae in total, 71 of which 
were observed with the CFHTLS. Their redshift range is be- 
tween 0.25 and 1. The remaining 44 objects are nearby super- 
novae, with 0.015 < z < 0.13. In the next section ( 12.2. U we 
describe the ingredients of the standard likelihood analysis, ig- 
noring systematics, which is similar to A06 with slight modifi- 
cations. Section l2.2.2l introduces the photometric calibration pa- 
rameters and defines the updated likelihood function including 
systematics. 





Fig. 3. Left panel: Id marginal for cr 8 for Q m = 0.25, the verti- 
cal lines are as in Fig. [2] Right panel: Correlation between co and 
<x 8 (68% and 95% confidence levels) for Q m = 0.25. The verti- 
cal lines indicate mean (dashed) and 68% confidence intervals 
(dotted). 



with unity but with large error margins, Co = 1.1 +0.6 (68%) (see 
also Fig. [2] right panel). At 95% confidence nearly the whole 
range of cq is permitted. Ignoring the high-z calibration error 
leads to smaller error bars for other parameters, in particular cr 8 
(left panel of Fig. For a fixed £2 m = 0.25 we respectively 
get cr 8 = 0.78^ and erg = 0.77 ± 0.05 with and without 
modelling the systematic error. Because we allow the calibra- 
tion to be smaller and greater than unity there is no significant 
bias on cr 8 . If we restrict cq to values smaller or equal 1, thus us- 
ing the information about the underestimation of the shear signal 
as prior, we obtain a best-fit <x 8 of 0.83 (Q. m = 0.25). Thus, cr 8 
might be underestimated by 8% which is about one sigma, pro- 
vided our simple model is correct. From Fig. [3] (left panel) we 
see the potential shear measurement systematics shifts the like- 
lihood and introduces a tail for high <x 8 . 

We repeat the above analysis with zo = 0.8 and 1.2 and 
find only small shifts A<x 8 for the power-spectrum normalisa- 
tion. Writing A<x 8 = A (zo - 1), we find A = -0.03 if co is varied 
between and 2, and A = -0.09 for cq < 1, corresponding in 
both cases to a 2% change for cr 8 in the considered range for zo- 

Interestingly, we see no correlation between cq and any other 
individual parameter when the whole posterior is considered. 
However, cq is of course correlated with the combination of Q m 
and cr 8 which determines the shear amplitude. When fixing D. m 
we see a strong correlation with cr 8 , as expected (Fig. [3] right 
panel). 

To summarise the results on weak lensing systematics, we 
note that this study is far from complete. Other tests focusing on 



2.2.1. Standard likelihood analysis 

We use the results from the SNIa light-curve fits of A06. For 
each supernova the fit yields the rest-frame B-band magnitude 
m* B , the shape or stretch parameter s, and the colour c. These 
quantities are combined in the estimator of the distance modulus 
for the i th object as follows: 



Mbj = m* Bi -M + a(sj - 1) -/fc,. 



(6) 



The universal absolute SNIa magnitude is M; a and f3 are uni- 
versal linear response factors to stretch and colour, respectively. 

The error needed for the likelihood includes: (1) The mea- 
sured uncertainties of the light-curve parameters m* B ., s,-, c, and 
their correlations. The corresponding error contribution to fi B j 
can be written as cr 2 (p B j) = 0^W 2 2 with the parameter vec- 
tor 02 = {l,a,/3) and the covariance W2 of the data vec- 
tor (m* B ., St, Ci). (2) An additional redshift uncertainty cr pv , = 
5/ In 10 ■ v p /(cz,) due to peculiar velocities corresponding to 
v p = 300 km s _1 . (3) An intrinsic scatter in absolute magnitude 
of cr int = 0.13. 

Those three errors are added in quadrature to the log- 
likelihood expression, which is 



xlniP) 



^,,(p)-51og 10 (^)] 2 



(7) 



The linear dependence on the Hubble constant of the luminosity 
distance c/l is taken out as an additional summand in eq. (0 and 
integrated into the absolute luminosity. From now on we will use 
the parameter M = M - 5 log 10 hjQ. 

The error term from the light-curve fit, cr 2 {p Bi ), depends on 
the stretch and colour responsitivies a and fj. A06 kept those 
parameters fixed during the ^-minimisation and updated them 
iteratively in consecutive minimisation runs. This leads to a bias 
in those parameters of 7% - 15%. Instead, we include the depen- 
dence of a and fj in the denominator in eq. (0 to obtain unbiased 
best-fit results for those parameters which are therefore larger 
than those cited in A06, see Table |4] 

We do not take into account the term involving the covari- 
ance determinant in the log-likelihood. Even though in this case 
this term is not a constant because of the parameter-dependent 
variance, we verified that the effect on cosmological parameters 
is very small. Moreover, it is in accordance to most o ther SNIa 
analyses (e.g. lAstier e t al. 20061 iKowalski et al.ll2008l) . 
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2.2.2. Extended analysis with systematics 



2.2.4. Bias due to systematics 



We extend the analysis of A06 by including the response of the 
distance modulus to a photometric zero-point shift in each of the 
seven filters (g,r,i,z,U,B,V) and in the Vega (B - R) colour 
in the Landolt system. The reference photometric zero-points 
are those used in A06, which were estimated without using the 
Hubble diagram. If 6\ denotes the vector of those eight zero- 
point magnitude shifts, the linear response kj a is the change in 
distance modulus for the 2 th supernova for a small change of the 
a th zero-point, 



1 A ^.' A n 

k <« = T7T~ ; tt = 0...7. 
A6i a 



(8) 



By changing the zero-points and redoing the light-curve fit, we 
obtain the values of k\ a for each supernova. We include this 
change to the distance modulus linearly in the model, which re- 
sults in an additional term k) ■ 6\ in the likelihood. 

We assume that the zero-point magnitude shift parameters 
6\ are uncorrected variables, since the data in different optical 
bands have been reduced independently. The noise is the result 
from a number of different source and can therefore be well ap- 
proximated to be Gaussian. By definition they mean of 6\ is zero. 
The rms is taken to be 0'. n 01 for all filters except for z (fin) where 
we assume the rms to be 0™03. The numerical values of those 
uncertainties are taken from A06 (Sect. 4.1) and correspond to 
the limits of reproducibility of the photometric calibration. The 
corresponding (diagonal) covariance matrix is W\ . 

This prior information is multiplied to the SNIa-likelihood 
in the form of a multi-variate Gaussian likelihood. With the ad- 
ditional term, the log-likelihood expression corresponding to the 
extended analysis is then 



z 



[^, < .( f> ) + t ).g 1 -51og 10 (^)] 2 

0~ 2 (lUBi) + V 2 ■ + O- 2 t 



(9) 



Now, the parameter vector p contains the zero-point parameters 
0\. The corresponding optical bands are indicated in Table [2] 



Table 2. List of SNIa systematic parameters and their symbols. 
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a 




Light-curve parameters 
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#23 




Ag 






Ar 


0u 




Ai 


012 




Az 


0,3 


Zero-point shifts 


AU 


014 




AB 


015 




AV 


016 




MB - fl) Vcg a 


017 





2.2.3. Systematics for SNIa 

Unlike systematics in weak lensing shape measurements which 
are difficult to model, the observation-related systematics for SN 
are more easily parametrised. In the next section we take into 
account errors in the estimated distance modulus due to uncer- 
tainties in the photometric calibration. The results of the MCMC 
analysis for the SNIa internal parameters are given in Table |4] 



We compare the case of ignoring the systematic errors, using 
eq. (0 as the log-likelihood, with the case of fully taking into 
account the systematics according to the log-likelihood (0. As 
can be seen in Fig. [4] ignoring the zero-point errors leads to 
an asymmetric decrease of the error bars. The constraints get 
tighter, mainly along the direction of constant luminosity dis- 
tance, which is the parameter-degeneracy direction. The error 
bars decrease by about 10%, see Table[3] The bias on parameter 
means is small, between 10% and 20% of the statistical uncer- 
tainty. The bias on the intrinsic SNIa parameters (M, a,B) is even 
smaller, not more than a few percent of the statistical uncertainty. 
For a fixed Q m = 0.25, the absolute biases on both the eos pa- 
rameter (for Ode = 1 - fi m ) and the dark-energy density (for 
w = - 1 ) are smaller than for the marginalised case, but remain 
to be about a tenth of the statistical error. 

As mentioned in Sect. 12.21 we obtain unbiased best-fit values 
for the stretch and colour response parameters, a and B, respec- 
tively (see Table|4]i. These differ by about 15% from the (biased) 
values given in A06. The absolute magnitude M is consistent 
with A06. 

The parameters describing the zero-point shifts (Oio to On) 
are all consistent with zero (Table 0J. Except for 613 = Az they 
have zero mean and rms of about 0™01. The influence of the 
derivative (eq. [8) on the second term of the likelihood (eq. |9]l 
is small in comparison with the first term. Nevertheless, corre- 
lations with other parameters are introduced as is discussed in 
the next section. The mean of the z-band zero-point shifts On is 
negative (although not significantly so) and its variance smaller 
than the expected value of 0'. n 03. It is also this parameter which 
shows the highest correlation with cosmological parameters. 

2.2.5. Correlation between systematic errors and 
cosmological parameters 

In most cases the zero-point shifts 6\ are uncorrected with other 
parameters. However, some pairs of one or more nuisance pa- 
rameters show correlations, most notably, the absolute SNIa 
magnitude M is correlated both with the B- and V-band uncer- 
tainty AB = 0i5 and AV = Ok,, respectively, see Fig. [5] This 
correlation can be explained by looking at the rest-frame colour 
parameter, c - (B — V)Bmnx + 0.057, which contains the B - V 
colour excess at the time of the B-band maximum. Redefining 



Table 3. Parameter means and 68%-confidence intervals when, 
respectively, ignoring and taking into account the systematic 
errors in the form of photometric zero-point fluctuations, see 
Sect. rZ24l 
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wCDM, 




0.2 0.4 0.6 



Fig. 4. Confidence contours (68%, 
95.5%, 99.7%) for full treatment of 
errors due to photometry zero-points 
(blue solid lines) and ignoring those 
errors (green dashed curves). These 
two cases correspond to the two cases 
in Table [3] The left panel corresponds 
to w = -1, the right panel is for a flat 
vvCDM Universe. The constraints are 
obtained using SNIa alone. 



Table 4. Mean and 68% errors for SNIa internal parameters, for 
a flat wCDM model. The values of On, i = . . .7 are in units of 
0.01 magnitudes. 



Parameter 


Best-fit-value 


-M 


1 Q 3T7+0.036 
ly - J -"-0.041 


a 


1.62!^ 


~P 




0io = Ag 


-0 06 +0 " 


On = Ar 


"• lu -0.96 


012 = a; 




013 = Az 


-1 +25 

1 ' -2.5 


14 = At/ 




0,5 = AS 


-0 04 +0 " 


016 = AV 


-0 06+ - 98 

"•""-0.97 


e„ = A(B- i/) Vega 


-0 05 +0 " 



the B- or V-band zero-points causes a systematic change in c 
which is compensated by a corresponding change in M. The 
slope of the measured M - AV-correlation is about 1.5. The cor- 
relation between M and AB = 9^ is much weaker and has a 
smaller slope of -0.3 to -0.4. The correlation between M and the 
difference AB - AV is therefore on the order of the best-fit value 
of B — 1.80, which is expected from eq. [6] The difference in 
correlation can be explained by the fact that the observed mag- 
nitudes of each object translate into different rest-frame bands 
for the light-curve fit depending on the object's redshift. 

Further, the z-band zero-point offset (On) shows a cor- 
relation with both £2 m and w (Fig. |6j. We find Aw = 
0.057 A0i 3 /(O.O3 mag) and AQ m = -0.0159 A0 13 /(O.O3 mag). 
The latter value can be compared to A06 (Table 5) who varied 
the zero-points to infer the influence on cosmological parame- 
ters; our value is of the same order of magnitude as the one from 
A06 but has opposite sign. For a fixed £2 m = 0.25 we obtain 
Aw = -0.0135 A0i 3 /(O.O3 mag), about half of the value cited in 
A06 (Table 5). 



CO 

< 





19.1 19.2 19.3 19.4 19.5 19.6 

-M 



19.1 19.2 19.3 19.4 19.5 19.6 

-M 



Fig. 5. Correlations between the universal absolute SNIa mag- 
nitude M and the fi-band (#15, left panel) and V-band (0i6, right 
panel), respectively. The zero-point shift parameters 6u are given 
in units of 0.01 magnitudes. 





< 

- 11 

- CO 



0.2 0.4 0.6 0.8 1.0 1.2 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

n m w 

Fig. 6. Correlations between the z-band zero-point offset (613) 
and D. m (left panel) respectively w (right panel). 



2.2.6. Further sources of systematics 

A06 discussed further sources of systematic uncertainties which 
can affect cosmological results from SNIa. Some of those un- 
certainties could be parametrised and invoked in a joint analysis 
with cosmology as presented in this paper. We state a few more 
examples here and leave a thorough analysis for future work. 

For example, the difference between measured rest-frame U - 
band magnitude and the one 'predicted' from the light-curve fit 
using only two bands can be included in the likelihood function. 
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In addition, the Malmquist bias can be modelled as a func- 
tion of redshift. The inclusion of this bias will be necessary when 
redshift-dependent dark-energy models are to be tested. This re- 
quires accurate knowledge of the sample and of selection effects. 

Furthermore, weak gravitational lensing of distant SNIa can 
be modelled using the dark-matter po wer spectrum, either fr om 
(non-linear) theoretical prescriptions (|M enard & Dalai 20051) or 
numerical simulations (Martel & Premadi 2008). 



2.3. WMAP5 CMB anisotropies 

To calculate CMB temperature and polarisation power- and 
cross-spectra we u se the publicly available package CAMf{] 
dLewis et al . 2000). The likelihood is evaluated using the public 
WMAP5 code! 8 ! dDunklev et alj|2009h . We include all 'standard' 
components which are the same for three- and five-year. Those 
are the low-^ TT and TE/EE/BB spectra, the high-/" TT and TE 
spectra, and the point-source TT correction. The low-f (I < 32 
for TT, I < 23 for polarisation) likelihoods are calculate d us- 
ing Gibbs sampling dPage et al.ll2007HDunklev et al.ll2009l). The 
high-l sampling uses pseudo-Q according to iHinshaw et al.l 
J2007I) . 

In contrast to the published WMAP5 results, we do not 
include corrections due to SZ. The SZ amplitude is uncon- 
strain ed by WMAP5 and is not degenerate with other parame- 
ters (iDunklev et al .1120091) . A recent analysis found no biases for 
WMAP5 whe n ignoring the contri bution of the thermal SZ effect 
from clusters dTaburet et al.ll2009l) . 



3. Parameter analysis and cosmological constraints 

3.1. Cosmological model and parameters 

We assume a flat dark-energy cold dark matter cosmology 
(wCDM) with the parameter vector (Q m , Q[,, t, w, n s , /z, <x 8 ). 
Clustering of dark energy is not taken into account. For CMB, 
instead of the normalisation defined at a scale of 8 Mpc/h, erg, 
the parameter which is sampled is A^, the curvature perturba- 
tions amplitude at the pi vot scale kg = 0.00 2 Mpc . As stated 
in Sect. [2J] we use the ISmith et all d2003l) fitting formula to 
model the non-linear power spectrum needed for weak lensing. 
We do not include a massive neutrino component to the mass- 
energy tensor and assume N s g = 3.04 as the effective num- 
ber of massless neutrinos dMangano et al.l 12002 ). which is the 
preferr ed value for WMAP5 lDunklev et al.ll2009l) . lTereno et alj 
(2008) obtained constraints on the neutrino mass using var- 
ious probes inclu ding CFH TLS-Wide w eak lensing, see also 
iGong et al.l d2008l) and llchiki etail d2009l) . Additional parame- 
ters are (a,b,c,co) from lensing and (M, a,/3, #10, ■ ■ ■ Qyj) from 
SN la. These parameters are described in the respective subsec- 
tions of Sect. |2j including systematics we sample a total of 22 
parameters. 



ITT 



3.2. Sampling the parameter space with MCMC 

We use an adapt i ve Metropolis- Hastings algorithm 
dMetropolis et all 119531; lHastingsl 1 1970b to generate Monte 
Carlo Markov Chains (MCMC) as sample of the posterior. 
The Fisher matrix evaluated at the maximum-likelihood (ML) 
parameter serves as initial multi-variate Gaussian proposal. 
The ML point is estimated by a conjugate-gradient search 



Lensing u/CDM - 
Lensing ACDM 
CMB wCBM 
CMB ACDM 



CO 

D 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Fig. 7. 68% and 95% confidence levels for lensing only, assum- 
ing wCDM (solid red lines), ACDM (dashed red), and CMB 
only, for wCDM (solid blue), ACDM (dashed blue curves). 



7 http://camb.info 

8 http://lambda.gsfc.nasa.gov 



dPress et alJ Il992l) . For lensing and SNIa the typical chain 
length is 100000 with an acceptance rate of about 15%-20%. 
Every 1 000 steps the chain covariance is updated and replaces 
the previous proposal. The proposal variance is multiplied by 
a factor 2.4/ V»dim which is optimal for a Gaussi an posterior 
dHanson & CunningharrJl998HDunklev et al.ll2005l) . 

For WMAP5 on the other hand, we choose a different strat- 
egy. With acceptance rates as stated above we found the poste- 
rior not well sampled and as a consequence the parameter un- 
certainties largely underestimated. We refrain from updating the 
proposal but kept the initial guess, which is larger than the Fisher 
matrix, for the whole duration of the chain. This choice results in 
small acceptance rates of 3% to 10% but guarantees an unbiased 
sampling of the posterior. 

The experiments are combined by multiplying the corre- 
sponding likelihoods or, equivalently, by summing up the log- 
likelihoods. The confidence intervals are obtained by creating 
histograms of the parameter vectors of the chain and estimating 
regions with 68%, 95% and 99.7% of the density. 

In particular for WMAP5 the calculation of the converged 
chain is very time-consuming and takes a few days to a week 
on a fast multi-core machine. To overcome this and other prob- 
lems related to MCMC, such as the difficulty to assess conver- 
gence and the high correlation of a Ma rkov chain, a new method 
called Population MonteCarlo (PMC, ICappe et al.ll2004 120071) 
has been developed which is based on adaptive importance sam- 
pling. In a companion paper, we test this method with simul a- 
tions and apply it to cosmology posteriors (Wra ith et al .1 12009). 

3.3. Combined constraints from lensing and CMB 

The main contribution from cosmic shear (without tomographic 
information), in combination with CMB anisot ropy measure- 
ments, are constraints on the normalisation erg dContaldi et al.l 
120031: iTereno et al.ll2005l) . For a ACDM model, the degener- 
acy direction between Q. m and <x 8 is orthogonal between lensing 
and CMB. Lensing is sensitive to the amount of structure and 
therefore, a high normalisation has to be countered by a smaller 
amount of the total matter. On the other hand, an overall increase 
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Fig. 8. 68% and 95% confidence levels for the three parame- 
ters which are affected most by systematics (O m , w, era). Solid 
(dashed) contours correspond to the case of ignoring (including) 
systematics. 



of the CMB angular power spectrum by a higher normalisation 
has to be compensated by more matter to dampen the peaks. 

This picture changes if the dark energy equation-of-state is 
a free parameter. Whereas the lensing contours for O m and <x 8 
do not broaden much, the CMB constraints increase dramati- 
cally. Most importantly, they increase along the lensing-direction 
of degeneracy, see Fig. Q The result is that adding lensing to 
CMB data will not improve the constraints in the n m -cr 8 plane 
as much as it would do in the framework of a cosmological con- 
stant model. 

Lensing without tomography is not very sensitive to w, 
thus leaving it a free parameter has a minor influence on the 
other parameters. For CMB however, strong correlations enter 
through the angular diameter distance which CMB measures 
only at z = 1100. Those degeneracies (see Fig. |9]l can be bro- 
ken by adding an additional distance mea surement at low red- 
shift as a lever, e.g. using BAO or SNIa. (jFrieman et al J 120031: 
iBlake & Glazebrookll2003l: ISeo & Eisensteinll2003l) ~ 

3.4. Combined constraints from lensing, SNIa and CMB 

We will now discuss the main results of this paper. Joint con- 
straints using lensing, SNIa and CMB are compared for the two 
cases with and without taking systematics into account. In the 
former, both lensing and supernovae systematics are included. 
The results are given in Table [5] Fig. [8] and Fig. [9] We show 
the best-fit angular and 3d power spectra in Fig. [10] A CAMB 
parameter file with our best-fit parameter values is available for 
downloacQ. 

With the current data, SNIa is more efficient than lensing 
in helping decrease uncertainties. CMB+SNIa gives nearly as 
tight constraints as CMB+SNIa+lensing for most parameters. 
The reason is that SNIa data show different degeneracy direc- 
tions, in particular for the pair Q m -w (see right panel of Fig. @). 



http://www2.iap.fr/users/kilbinge/params.ini 



This helps to pin down w and thus, the main degeneracy for 
CMB is largely lifted. As stated in the previous section, this is 
not the case for lensing — lensing without tomography cannot 
constrain w. The consequence is that even parameters to which 
SNIa is not sensitive, e.g. Qb, are very accurately determined for 
the combination CMB+SNIa. 

Nevertheless, lensing improves constraints from WMAP5 
substantially. Some CMB-related near-degeneracies which arise 
in the vvCDM model are partially lifted in combination with lens- 
ing. In particular for (£2 m , cr 8 ), there is a large gain when lensing 
data is added. 

The effect of systematics on the parameters means and errors 
can be assessed by comparing the last two columns of Table [5] 
The shift of the best-fit values is less than 15% of the statistical 
error in all cases. Including systematics in the analysis increases 
the error bars by 10%-35%, where £2 m , w and cr 8 are affected 
most. Varying the lensing redshift for the bias zo from 1.0 to 0.8 
and 1 .2 changes the results by less then one percent. 

4. Discussion 

In this paper we combine three different cosmological probes to 
test a possible deviation from a cosmological constant: (1) Weak 
gravitational lensing as a probe of structure formation and ge- 
ometry in the redshift range of about 0.2 to 0.8. (2) Supernovae 
la as standard(isable) candles up to redshift unity. (3) CMB 
anisotropies supplying a wealth of information of the recombi- 
nation era (z ~ 1100) and, to a lesser extend, of the Universe 
up to the re-ionisafion epoch (z ~ 6 - 15). The data sets are (1) 
CFHTLS-Wide T0003 (see F08), (2) SNLS first-year (A06) and 
(3) WMAP five-year (H09), respectively. 

We test models in the context of a flat CDM cosmology 
with a dark- or vacuum-energy component with free but con- 
stant equation-of-state parameter w = pipe 1 . This corresponds 
to the simplest extension of the 'vanilla' ACDM model which 
goes beyond a cosmological constant. 

The joint constraints including the modelled systematics in 
the data sets are w = -1.02±°^+°-}* (68% and 95% confidence, 
respectively). Without taking the systematics into account, the 
resultisw = -1.01 ±0.06 + 0. 12, representing 25% smaller error 
bars. With the current data there is no evidence for a dynamical 
dark-energy component not being the cosmological constant. 

Two potential sources of bias in the third-year CFHTLS- 
Wide lensing data are scrutinised. One, the measured variations 
of the shear signal (aperture-mass dispersion (M% p )) between 
MegaCam pointings are compared to /V-body simulations. We 
estimate the measured fluctuations to be higher than expected by 
not more than about 5% to 15% on scales below 30 arc min- 
utes. Whereas this might be a hint of systematics in the data, it 
is not straightforward to model and to assess its effect on cos- 
mology. The second issue are systematics in the shape measure- 
ments which seem to lead to an underestimation of the lens- 
ing signal at high redshift. We devised a very simple model of 
this potential systematics by multiplying the lensing efficiency 
above z = 1 with a constant cq > 0, to mimic the effect of 
a decreased measured shear. Marginalising over cosmological 
parameters, using weak lensing alone, yields co = 1.1 ± 0.6. 
All probes combined do not constrain cq much better, we find 
Co — 1.1 + 0.5. Restricting ourselves to cq < 1, implying an un- 
derestimation of the lensing signal, increases cr g by about 8% for 
a fixed £2 m = 0.25, which is roughly equal to the statistical error. 
Therefore, in the framework of this simple model, the value of 
crg might be biased by 8% if this effect is ignored (as it has been 
in F08). 
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Table 5. CMB, lensing and SNIa in various combinations. The mean and 68% marginals are given. For the first four cases system- 
atics are ignored, the last column includes all systematics, from both lensing and supernovae (see Sect. |2.2.4| i. 



Parameter 


CMB 


CMB+Lens 


CMB+SN 


CMB+Lens+SN 


CMB+Lens+SN+sys 


( \ 


n n/ic+0.020 

11045 -0.016 


fl C\A 1 +0.016 


n rvin+0.0028 
0.043 j 0()026 


n nzn? +0 - 0026 

"■™ -0.M23 


n n/i9c -i- n aa9q 


o 

ai m 


A 9^9+0.099 

o.zoz _ 0O93 


fl 9/19+0.092 
U - /42 -0.048 


a 9^7+0.025 
' -0.023 


a 9CO+0.018 
-0.016 


n TCI +0.023 
-0.01! 


T 


0.087 ± 0.016 


0.086+°™ 


0-088!°,™ 


0.088+°;™ 


0.088 ± 0.017 


w 


-1.08_*» 


-1 09+ 024 

-0.22 


-1.025.°;™ 


-i-oio+^ 


-1 021+ a ° 79 

-0.081 


n s 




0.961*™ 


0.962 ±0.015 


963+ 015 

"■-'"-'-0.014 


963+ 014 




9 4^+0.13 


2.418!^ 


9 4^+0.12 


9 4 |4+o.098 


2.41 ±0.11 


h 


0.74+°!* 


a 7^4+0.096 


0.719+°;°| 


720+ 023 

-0.021 


723+ - 027 

-0.025 


erg 


0.82+°;^ 


8 19+° 061 

"■ 01 -0.069 


0.807+°;™ 


A 7QC+0.030 
U,/:,J -0.0Z7 


798+ - 037 

,;,o -0.044 



Combining all probes we assess the influence of systematics 
on the cosmological results. The contribution of systematics to 
the total error budget ranges between 5% and 40% (Table[5]l. The 
parameters which are affected most by systematics are Q m , a% 
and w for which this contribution is greater than 20%. Those are 
the parameters for which weak lensing adds significantly to the 
joint constraints. 

There are indications of more unaccounted systematics in 
the lensing data. Discussions about the origin of those system- 
atics and ways to remove them are addressed elsewhere (van 
Waerbeke et al. in prep.). Our findings strengthen the confi- 
dence in cosmological results from cosmic shear, presented in 
this work and earlier, using the same data, in F08. 

Our constraints are slightly tighter than ISeliak et"all d2005l) 
who obtained w = -0.99^^ using WMAP1, SNIa and, 
as probes of structure formation, SDSS galaxy correlations 
and Lyg forest clustering . Similar constraints were quoted in 
IWang & Mukheried (l2006h . w = -1.00+^|+o.J6 from WMAP3, 
SNL S and SPSS. Our re sults are comparable to the ones stated 
bv lKomatsu etafl (T2009). -0.11 < l+w < 0.14(95%) stemming 
from WMAP5, BAO and SNIa. Consi s tency is also achieved 
with the constraints from iMantz et"ail d2008l) who combined 
WMAP5 with the X-ray cluster mass function, cluster baryon 
fraction and SNIa to get w = -1.02 ± 0.06 (68.3%). 

Earlier results including weak gravitational lensing used less 
wide and/or more sh allow data. Using th e first-year release of 
the CFHTLS (TOOOl. lHoekstra et alj2006l) . and with rather tight 
priors on other cosmologica l parameters, an upper bound of 
w < -0.8 (68%) was derived. ljarvis et al (2006) combined weak 
lensing from the 75 deg 2 CTIO survey with SNIa and CMB and 
obtained w = -0.89+° ^ (95%). The error estimate for the latter 
includes statistical uncertainty, PSF systematics and shape mea- 
surement calibration biases. 

This paper puts particular emphasis on the treatment of sys- 
tematics in the data. Known observation-related systematics are 
parametrised if possible and included in the analysis together 
with cosmological parameters. This allows us to directly quan- 
tify the influence of systematics on cosmology and to find pos- 
sible correlations. For SNIa, we calculate the response of the 
distance modulus to fluctuations of the photometry zero-points. 
Ignoring these error sources leads to an underestimation of the 
parameter errors by 10%. The parameter values itself are biased 
by a fraction of 10% to 20% of the statistical uncertainty. 



5. Outlook 

For future, high-precision experiments it will be of great impor- 
tance to understand and control systematic effects. In particular 
for exotic models, such as time-varying dark-energy or modifi- 
cations of GR, more and more subtle influences of systematics 
mimicking a signal have to be excluded. This calls for combined 
analyses of cosmological and systematics parameters. 

Weak lensing serves as an important and independent probe 
of cr 8 . The measurement of the lensing skewness will allow us to 
lift degen eracies between parameters, in particular between Q. m 
and cr 8 ( Kilbinger & Schneider! 120051) . Power and bi-spectrum 
tomography can be used to constrain the time-e volution of the 
dark-energy eos parameter ( Taka da & Jainl [2004b . In addition, 
as a complementary and nearly independent probe of the non- 
Gaussian, high-density regime of the large-scale structure, weak 
lensing cl uster counts can imp rove dark-energy parameter con- 
straints d^kadi]&Bridli[2003). 

For lensing, systematic effects which have to be taken into 
account for future analyses include measurement errors (galaxy 
shapes, PSF correction, photo-z's), astrophysical sources (intrin- 
sic alignments, source clustering) and theoretical uncertainties 
(non-linear and baryonic physics) . The modelling of all those 
effects, as suggested recently by iBernsteinl d2008l) . leads to a 
huge number of parameters, on the order of several dozens to 
hundreds, depending on the number of redshift bins. In case of 
simultaneous determination of cosmic shear and intrinsic align- 
ments a high number of redshift bi ns is required (Bri dle & Kind 
l2007t iJoachimi & Schneiderll2008l) . 

Already with the present data, more general models of dark 
energy can be constrained. A wide variety of s uch models are 
tested with recent SNIa data (Ru bin et alJl2008l) . Because dark 
energy and modified gravity can influence distances in a non- 
distinguishable way, it is important to include measures of the 
growth of structure. A survey like the CFHTLS-Wide reaches 
out into the linear regime and is less prone to small-scale un- 
certainties due to baryonic and non-linear physics, which makes 
such a survey an excellent probe of the present acceleration of 
the Universe. 
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